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ABSTRACT 

We explore the use of Integral Field Spectroscopy (IFS) for observing extrasolar planet transits. 
Although this technique should find its full potential in space observations (e.g. JWST, TPF), we 
have tested its basics with ground based time series observations of HD209458b obtained with WHT 
+ INTEGRAL during a transit in August 17/18, 2004. For this analysis we have used 5550 spectra 
(from a potential of ^30000), obtained in 150 exposures during a period of more than 7 hours. We 
have found that IFS offers three fundamental advantages with respect to previously used methods 
(based on imaging or standard slit spectroscopy). First, it improves the effective S/N in photon 
limited observations by distributing the light coming from the star into the two dimensions of the 
detector. Second, this type of IFS data allows to autocalibrate instrumental and background effects. 
Third, since the star image characteristics (i.e. seeing, spatial shifts, etc.) as well as its photometric 
properties are extracted from the same data-cube, it is possible to decorrelate photometric instabilities 
induced by PSF variations. 

These data have also allowed us to explore the accuracy limits of ground based relative spectropho- 
tometry. This was done using a photometric index that probes the NaD lines, for which we obtained 
a nominal 1-a error of ~1.0 • 10~ 4 . This result, based on observations of only one transit, indicates 
that this type of ground observation can constrain the characterization of the transmission spectrum 
of extrasolar planets, especially if they cover multiple transits under good weather conditions. The 
present observations are compatible with no extra NaD depression during the transit. Although this 
result seems to be inconsistent with the recently reported HST - STIS findings we point out its limited 
statistical meaning: the results disagree within 1-a, but agree within 2-a. 

This method requires careful and specific reductions, and details about this process are given. We 
also give some recommendations to instrument developers in order to enhance the efficiency of the 
method. 

Subject headings: techniques: integral field spectroscopy, planetary systems, stars: individual 
(HD209458) 



1. INTRODUCTION 

The transit method (TM) has recently gained much 
attention among the various methods for the detection 
of extrasolar planets due to its potential for discovering 
earth-like planets. This method also offers the possi- 
bility for studying the existence of planetary satellites, 
rings, and atmospheric features (see, for instance, Deeg 
2002). All these possibilities have been clearly demon- 
strated with HD209458, the first planetary transit dis- 
covered (Charbonneau et al. 2000; Henry et al. 2000; 
Mazeh et al. 2000). The light curve of this object has 
allowed an unprecedented characterization of the stel- 
lar/planet system (e.g. planetary and stellar radii, stel- 
lar limb darkening, orbital inclination, period, etc), im- 
posing limits to the presence of planetary satellites and 
rings (Brown et al. 2001), and detecting features asso- 
ciated with the planet atmosphere (Charbonneau et al. 
2002; Vidal-Madjar et al. 2003, 2004). 
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The success of the TM largely relies on the (spectro) 
photometric accuracy for which the light curve can be 
obtained. For some science applications (e.g. detection 
of planetary atmospheric features) only accurate rela- 
tive spectrophotometry (i.e. bandl/(nearby)band2 and 
in/out transit) is required, reducing drastically system- 
atic correlated errors which may otherwise limit the S/N. 
Experiments with the HST (Gilliland, Goudfrooij, Kim- 
ble 1999; Brown et al. 2001; Gilliland and Arribas 2003) 
have demonstrated that accuracies of ~10 -5 or higher 
can be reached from space. 

Such a high accuracy is not always required. For in- 
stance in the wavelength domain 1190-1710 A, Lyman a, 
OI, and CII lines suffer during the transit a drop of ^15 
, 13 and 7.5 percent, respectively (Vidal-Madjar et al. 
2003, 2004). In optical wavelengths, the variations are 
indeed smaller. Charbonneau et al. (2002) have found a 
drop of 2.5 parts in 10000 (i.e. 0.025 percent) during the 
transit in the NaD resonant lines. This drop is smaller 
than expected from (poorly constrained) models, which 
also predict that other spectral characteristics in the op- 
tical (e.g. Rayleigh continuum) should be detected with 
lower accuracies (e.g. Brown 2001). The smaller than ex- 
pected NaD absorption detected by Charbonneau et al. 
(2002) may have resulted from a high cloud deck, a low 
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atomic Na abundance, or a combination of both. 

In this context, it is relevant to understand which are 
the accuracy limits that can be achieved from the ground. 
This should allow us to know if part of future studies can 
be performed taking advantage of the more abundant and 
flexible observing time available from ground-based ob- 
servatories, or if we need to rely entirely on space data. 
Note that, with the advent of large aperture ground tele- 
scopes, the Poisson limit may be boosted allowing access 
to more distant systems. This is particularly relevant 
in connection with follow up observations of the plane- 
tary transits to be discovered by missions such as Ke- 
pler (Borucki et al. 2004). Pioneering attempts to detect 
spectral features due to the planetary atmosphere from 
the ground in 51 Peg and HD209458 include those of 
Coustenis et al. (1998), Rauer et al. (2000), Bundy & 
Marcy (2000), Moutou et al. (2001), and Brown et al. 
(2002). More recently Deming et al. (2005) using NIR- 
SPEC on the Keck II telescope have reported upper lim- 
its to the CO absorption at a sensitivity level which sug- 
gests that a general masking mechanism is present in 
the planetary atmosphere of HD209458b. Also recently 
Narita et al. (2005) have studied the transmission spec- 
trum of this object on the basis of high resolution spec- 
troscopy with Subaru. Although the accuracy of these 
observations were unable to confirm or contradict the re- 
sult of Charbonneau et al., they reported upper limits on 
absorption due to several optical transitions. All these 
ground based observations have gained importance since 
STIS-HST is no longer operational. 

In this paper we propose to take advantage of the com- 
prehensive collecting nature of IFS to perform this type 
of observation. The observational method, whose advan- 
tages are described in Section 2, is tested with ground 
based data. The observations are described in Section 3. 
Details about the reduction process for this particular 
application are described in Section 4. The results are 
shown and discussed in Section 5, while Section 6 sum- 
marizes the main conclusions. 

2. PROPOSED METHOD: THE POTENTIAL OF IFS IN 
OBSERVATIONS OF PLANET TRANSITS 

Gilliland, Goudfrooij, and Kimble (1999) proposed to 
increase the photometric accuracy of very high S/N ob- 
servations by using the spectroscopic mode of STIS. By 
distributing the star light along the dispersion direction 
it is possible to increase the total number of photons col- 
lected before the detector reaches the saturation limit. 
This simple idea enhances the duty cycle of the instru- 
ment (i.e. total number of photons by unit of available 
time, including overheads), which translates to an in- 
crease in the S/N when the photon noise is the major 
source of noise. The larger detector footprint of the stel- 
lar light has a negligible consequence in terms of adding 
noise in a very high S/N regime. Actually it has the 
beneficial effect of averaging out flat field and sensitivity 
residuals. Another practical advantage is the possibility 
to define a posteriori ad hoc 'filters' well suited for the 
science applications. 

Here we propose to extend this idea by means of in- 
tegral field spectroscopy (IFS). Thanks to the 2D — > ID 
conversion done by a fiber bundle (or by an image slicer), 
the star light is not only spread along the dispersion di- 
rection but also across dispersion. This technique allows 



us to use a larger area of the detector, boosting the limit 
imposed by the photon-noise per exposure. Note that to 
optimize the use of the detector, the image of the star at 
the focal plane should be well oversampled by the Inte- 
gral Field Unit (IFU). Since for these transit observations 
we are interested in the total number of photons recorded 
in a limited period of time (i.e. duration of the transit), 
the image at the input of IFU could be defocused. Note 
that defocussing the star image on the IFU does not mod- 
ify the spectral resolution of the observations, which is 
another practical advantage of the IFS method with re- 
spect to previously used spectroscopic methods. 

The ability to distribute the light over the whole 
detector is particularly relevant with the advent (and 
prospects) of very large telescopes. This can be easily 
understood with the following example. Typical imaging 
observations of HD209458b (V = 7.64) in a 100 A band 
with an 8m. telescope will saturate in ^0.1 sec. Consid- 
ering a typical readout time of 60 sec, and taking into 
account that the transit lasts ~3 hours, an effective in- 
tegration time of ~20 seconds on the source during the 
transit could be obtained. If the observations are now 
done with a standard slit spectrograph working at, say, 
1 A/pixel, saturation is roughly reached in ~10 sec, and 
the effective integration time becomes ~30 minutes. If 
now the observations are done with an IFS at the same 
spectral resolution and sampling the (defocused) star im- 
age with, say, 100 spatial elements on the IFU (spaxels, 
according to IFS community), the effective exposure time 
may be increased up to 2.8 hours (out of 3 available). 
Of course, real observations can be optimized by defin- 
ing faster readouts, modifying the spectral resolution, 
etc. Detectors like that used by ULTRACAM (Dhillon & 
Marsch 2002) may also reduce overheads. The numbers 
above are just illustrative of current standard detectors. 

Apart from the ability to gather photons, stability is 
crucial for TM observations. Of course, beyond some 
limits it is not possible to control the stability of the en- 
vironment and/or the observational set-up. While this 
is true for both space and ground observations, the un- 
stable nature of the Earth atmosphere makes very accu- 
rate ground observations particularly challenging. In any 
case, when uncontrolled instabilities may affect the pho- 
tometric accuracy of the system, a good record of them 
may give the possibility to decorrelate the photometric 
signal. This provides another important advantage of 
the IFS. Since both, the spectra from which the pho- 
tometric information is obtained and the images of the 
object (PSFs) are extracted from the same data cube, 
we can remove photometric variations induced by PSFs 
instabilities. Furthermore, in fiber-based IFS (and to 
a lesser degree a slice-based IFU) a shift of the image 
on the focal plane may be well tracked independently of 
a shift between the pseudo-slit and the detector. This 
gives full control of the PSF variations, as well as the 
optical and mechanical shifts during time series obser- 
vations. Fiber systems also provide a series of generic 
advantages such as, for instance, the azimuthal scram- 
bling of the light (which reduces errors in the baryecnter 
of the recorded fiber image), the possibility to use static 
spectrographs (which reduces flexure induced errors), etc 
(see more details in, for instance, Arribas & Mediavilla, 
2000). On the other hand, potential problematic aspects 
of fiber-based IFS systems are the modal noise due to 
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small changes in fiber position or illumination (e.g. Bau- 
drand and Walker 2001), and a possible temporal varia- 
tion of the fiber throughput with wavelength. However, 
we will see later that, at least for the present instru- 
mental configuration, our data do not suggest they are 
significant sources of noise. 

Another advantage of IFS for these type of observa- 
tions is a consequence of the large number of spectra 
collected simultaneously, which can be used to autocali- 
brate the data themselves from detector and background 
signatures. This will be shown in Section 5. 

Summarizing, IFS provides three main advantages for 
transiting planets observations: i) improves the ability 
to collect photons during the transit, enhancing the S/N, 
ii) it allows to autocalibrate the data, and iii) by tracking 
in an independent manner the instabilities produced at 
the focal plane and those due to the spectrograph, it is 
possible to remove noise correlated to PSF characteris- 
tics. 

3. OBSERVATIONS 

We performed the observations with the 4.2m William 
Herschel Telescope during 2004 August 17/18 using the 
INTEGRAL system (Arribas et al. 1998) connected to 
WYFFOS (Wide Field Fiber Optic Spectrograph; Bing- 
ham et al. 1994). We used the recently installed WYF- 
FOS long camera, equipped with two EEV-42-80 thinned 
and coated CCDs butted along their long axis to pro- 
vide a 4K x 4K pixel mosaic. The high speed read out 
mode (45 sec) was used, with gains of 2.3 e- (chipl-red) 
and 1.94 e- (chip2-blue) yielding a readout noise of 6 
e-. These detectors are considered linear for standard 
observations until ~52 K ADU (see below). Note that 
WYFFOS is mounted in one of the Nasmyth platforms 
of the WHT, which provides good mechanical stability 
(see below). 

The bundle of fibers SB1 (Standard Bundle 1), consist- 
ing of 205 hbers, each 0.45" in diameter, was used. This 
bundle is arranged such that 175 hbers cover a rectangu- 
lar area of 7.8" x 6.4" while 30 additional hbers, forming 
a ring of 90" in diameter concentric with the rectangle, 
measure simultaneously the sky background. The actual 
distribution of the hbers at the focal plane can be found 
in Figure 6 of Arribas et al. (1998) and also in Figure 1 
of del Burgo et al. (2000). 

The spectra were taken using the R316R line/mm grat- 
ing with an effective resolution of ~3 A and covering the 
-4000-10000 A spectral range. 

INTEGRAL does not have an independent focus mech- 
anism for the hber bundles, so the only possibility to 
spread the light over the input of the fiber bundle in a 
controlled way was to defocus the telescope. However we 
decided not to apply any major defocus to the telescope, 
because it would have affected the guiding quality. 

The observations on target consisted of 8 time series 
of 20 individual exposures of 120 sec each (the 8th series 
had to be suspended after 10 individual exposures were 
taken due to the detectable flux increase due to the 
twilight). After each time series we took a short flatfield 
exposure. Therefore, the total object data set consisted 
of 150 individual exposures (with a potential of —30000 
individual spectra though, of course, many had low S/N), 
extending for a period of more than 7 hours, from which 
5 hours were integrating on the target. The weather 



conditions were good, with a stable seeing of —0.8 arc- 
sec. The values of atmospheric extinction in the visible 
(r'-band) obtained by the Carlsberg Meridian Telescope 
(http://www.ast.cam.ac.uk/~dwe/SRF/camc_extinction.html 
also suggest stable and photometric conditions for the 
night, with a mean value for the extinction coefficient 
of 0.091 ± 0.011 magnitudes per airmass using data 
expanding a period of 4 hours. Non-photometric data 
were not detected this night. 

We located the object slightly off-center (i.e. around 
fiber number 111 in the standard code for this bundle) in 
order to minimize the effects of one broken fiber located 
in the opposite direction. We will refer to it as central 
fiber or fiber 111, and it is the closest to the peak of the 
image in most of the exposures. The 37 fibers around it 
covered most of the target flux. 

Apart from lamp-flats taken between sets of 20 expo- 
sures, bias, wavelength calibration, and sky flat images 
for calibration purposes were taken at the beginning and 
at the end during the twilight. The internal calibration il- 
lumination does not mimic the exit pupil of the telescope, 
producing slightly larger hber images. This should be ir- 
relevant provided the use of the these calibration images 
for the present application (see next section). 

4. REDUCTIONS 

The reductions of these IFS data have the goal of gen- 
erating both spectra (from which the relative spectropho- 
tometric index is derived) and images (from which PSF 
characteristics are obtained) for each of the 150 individ- 
ual exposures. Here we will briefly describe the steps 
followed during this process, commenting on the relative 
importance of each step for this particular application. 

4.1. Spectral reductions 

The present study will be focussed on the analysis of 
the NaD lines. As mentioned in Section 1 these lines have 
been studied with great accuracy with HST (Charbon- 
neau et al. 2002), which provide us with a good reference 
to compare the present observations. For this reason, 
and with the aim of simplifying and making the reduc- 
tion process more manageable, we trimmed our original 
images into a spectral region of ~200 A including the 
above mentioned resonant lines. 

4.1.1. Bias: 

The bias value was obtained from the overscan strips 
signal. For the red chip used here, the mean value in 
that region increased linearly over the night from 1018 to 
1028 ADUs. The rms of this relation is about < 1 ADUs, 
a very small value compared with the typical signal of the 
brightest spectrum (~40000 ADU). This indicates that 
the expected errors associated with an incorrect bias sub- 
traction should have a minimal effect on the relative pho- 
tometry of these very high signal to noise observations. 

4.1.2. Flat Field correction: 

By construction, in many IFS systems is not possi- 
ble to illuminate directly the detector with a uniform 
source (e.g. blank sky). This is due to the presence of 
the IFU (i.e. hber bundle, image sheer, microlens array) 
in the optical path. This prevents obtaining directly flats 
fields for correcting pixel to pixel sensitivity variations. 
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Instead flats usually refers to images obtained when illu- 
minating uniformly the input of the IFU, and sometimes 
are called flat-spectra. These are generally used to cor- 
rect both, differences in throughput among the different 
sampling elements (i.e. fibers), and sensitivity variations 
of extracted pixels. Note that an extracted pixel is gen- 
erally obtained by integrating several individual pixels 
across the dispersion direction, so they represent an av- 
eraged (weighted) behavior. Although this approach is 
generally enough for most applications, in some cases it 
may represent a limitation. We will show below that 
our spectra drift about 0.6 detector pixels during the 7.5 
hours of the observations, so different pixel sensitivities 
could emulate object flux variations. 

Although it is not possible to obtain direct detector- 
flats with INTEGRAL/WYFFOS, its detector is also 
used at the prime focus of the WHT for direct imag- 
ing. Fortunately this detector was used for several imag- 
ing programs during September 2004. So we collected 
the flats images recorded during this period (kindly pro- 
vided to us by D. Lcnnon) to create a very high S/N 
super flat. For the spectral range considered here we com- 
bined about 40 images (mostly in the V filter) and cor- 
rected them from low spatial frequency variations using 
the IRAF tasks imcombine and flat Id, respectively. Note 
that although relatively large fringing amplitudes have 
been measured at longer wavelengths (A > 7000 A), the 
wavelength range analyzed in this paper is free of that 
effect. This superflat has about 10 6 ADU ( 2.3 • 10 6 e") 
and, therefore, a nominal S/N per pixel of ^1500. 

4.1.3. Linearity correction: 

The photometric index at a given time is obtained after 
averaging the individual values obtained from the spec- 
tra recorded at that particular time. The fact that the 
relative intensity of these individual spectra may change 
during the night (due not only to guiding and seeing in- 
stabilities, but also to the unavoidable differential atmo- 
spheric diffraction), makes the spectral index sensitive to 
detector linearity deviations. 

The linearity of the present detectors was studied by 
Tulloch (2001). We have fit a polynomial function to 
his experimental data with signals <38 K ADU (rms 
of 0.14 percent). For signals >38 K ADU, the slope of 
the linearity curve changes sharply from its relatively 
smooth behavior at lower signals. This slope change is 
difficult to model in detail, so we decided not to correct 
for this linearity slope change (i.e. we apply a constant 
correction for values >38 K ADU). Linearity corrections 
at the lower signal end (i.e. <1000 ADU) are also difficult 
to model with high accuracy, but effects on the relative 
photometry are small. In Section 5, we will analyze and 
correct for the residuals of this linearity correction. 

4.1.4. Extraction: 

The detector readout noise has a minimal effect in 
these very high signal to noise observations, we there- 
fore decided to use relatively wide extraction windows 
of 14 pixels. This width guarantees the collection of 
most of the flux from a fiber/spectrum, but still avoids 
overlapping information from different spectra. In addi- 
tion, taking into account the relatively small drift of the 
spectra in cross dispersion over the night (~0.6 pixels 



in 7.5 hours), we decided to use a fixed extraction win- 
dow for all the spectra/exposures. (However, we also re- 
duced the data allowing for global shifts in the location 
of the extraction windows, but no appreciable changes 
were detected.) The functions defining the centers of the 
extraction windows were obtained by adjusting polyno- 
mials to the spectra of the image obtained by adding all 
the flat-lamp exposures (after previous removal of the 
cosmic rays). It was checked that these spectra were ac- 
tually located near the mean position during the night. 
These functions were 5th degree polynomials obtained 
after fitting 130 data points along the selected spectral 
range. Typical rms were of 0.015 pixels (no significant 
improvement was found when using higher degree poly- 
nomials). The extraction was done as a simple sum- 
mation (i.e. not with PSF shape and inverse variance 
weighting) in order to avoid undesired renormalization 
factors which could affect the flux. Although these fac- 
tors should have a very small effect in standard observa- 
tions, one could think that they may have some relevance 
for the current ultra-precise photometry. Therefore, we 
decide not to use this approach which, in addition, has 
very little advantage in a high S/N regime. In any case, 
we checked that no relevant differences were found when 
an optimal extraction was performed. 

4.1.5. Wavelength calibration: 

For the present analysis the wavelength calibration is 
important to define accurately the bands used for creat- 
ing the spectral index. Although the absolute conversion 
from pixels to wavelength coordinates does not need to 
be particularly accurate, it is important to correct any 
drift or change of the wavelength solution among fibers 
and, especially, during the night. 

For obtaining the absolute wavelength calibration for 
each fiber we used the wavelength calibration exposures 
obtained for such purpose. That was done using as ref- 
erence the brightest spectrum with the IRAF task cal- 
ibrate. Then, the individual wavelength solutions for 
each fiber were obtained with recalibrate. To detect any 
drift along the dispersion direction we cross-correlated 
the three brightest spectra in the images obtained along 
the night. These spectra are located in the region of 
the detector where most of the signal is extracted. We 
find a similar behavior in the three cases with a linear 
shift of ~0.3 pixels over the night at a constant rate of 
^0.04 pix/h, with a typical rms of 0.05 pixels. Therefore, 
to improve the accuracy of the drift associated with each 
exposure we combined all the spectra obtained in a given 
exposure, and cross-correlated these integrated spectra. 

Once the wavelength solution (and spectral drift) was 
known, we decided not to resample the individual spectra 
to a common spectral range and dispersion. This step 
would have implied an extra interpolation of the data. 
We found it preferable to use this information to define 
homogeneously the spectral bands of all the spectra (see 
Section 5.2). The wavelength solution was applied with 
the task dispcor of IRAF. 

4.1.6. Cosmic ray rejection: 

The cosmic ray rejection is an important step in the 
reductions, since a relatively small number of badly re- 
moved cosmic rays could degrade substantially the S/N 
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of the data. After trying several automatic rejection 
methods we decided that none of these methods was suf- 
ficiently reliable for our particular case. The reason for 
this is that, because the spectra are well separated on 
the detector, slight shifts of the image at the focal plane 
may produce relatively important changes in the flux at 
the detector, which may be interpreted as CR when dif- 
ferent exposures are compared. Therefore we decided to 
follow the secure method of removing the CR in an inter- 
active manner. In order to identify more easily cosmic 
rays, each individual spectrum was divided by a high 
S/N clean template, so even if a CR hits in the middle 
of the spectral feature (i.e. the NaD lines) it could be re- 
moved reliably. After the cosmic rays were identified and 
removed the spectra were multiplied by the template. A 
total of ~80 cosmic rays were identified and removed in 
the relevant spectral range but, of course, the full set of 
over 5000 spectra were inspected. 

4.1.7. Throughput correction: 

Flat-spectra were applied to correct for fiber sensitiv- 
ity. This was done using the blank sky observations dur- 
ing the twilight at the beginning and at the end of the 
night. Although the fiber-to-fiber sensitivity variations 
are large (rms ~20 percent) these are stable during the 
night, with typical differences of smaller than ~1 percent 
(i.e. the difference of two flats taken at the beginning 
and at the end of the night gave a rms of 1.1 percent 
for the 37 fibers used, which probably imposes an up- 
per limit to the actual fiber throughput stability). Note 
that global fiber sensitivity variations have no impact on 
the determination of spectral indices, which always probe 
relative flux variations (i.e. bandl/band2). A temporal 
variation in the fiber throughput with wavelength could 
be a limitation, but we don't have evidence for such a 
behavior. (As we will see in Section 5, changes in the 
throughput-wavelength function seem to be dominated 
by atmospheric more than for instrumental effects). 

Although the photometric index is not dependent on 
the correction of fiber throughputs, this affects the star 
image reconstruction. In fact errors in the fiber through- 
put corrections may produce incorrect PSF reconstruc- 
tion, introducing further uncertainties during the decor- 
relation process. The throughput correction may also af- 
fect the background subtraction. However, in these two 
cases the effects on the final results are very small (see 
below). 

4.1.8. Background subtraction: 

As described in Section 3, 30 fibers located in a ring 
of 90" in diameter around the star allowed us to obtain 
contemporaneous background information. The spectra 
of some of these fibers are located relatively close to fibers 
with very high signal, so they may be affected by cross- 
talk (or scattered light). Therefore, the averaged value 
of only 18 of these fibers was actually used to obtain the 
background information. 

4.2. Reconstructed images 

In order to obtain the star images, we created files 
with the X and Y positions of the fibers in the focal 
plane of the telescope and the integrated (background- 
subtracted) flux in the range 5829.200-5957.321 A used 



for the definition of the photometric index in Section 5 
(i.e. ft = h + h + fs)- With the help of the NAG rou- 
tine E01SAF, an interpolating 2-dimensional surface F(x, 
y) is generated. This routine guarantees that the con- 
structed surface is continuous and has continuous first 
derivatives. The interpolant F(x, y) was then evaluated 
regularly each 0.05" with the routine E01SBF on a grid 
of 80 x 80 pixels to create the maps presented here. Only 
37 fibers were considered to have enough signal for the 
definition of the photometric index (see Section 5.2), and 
these were also the ones used for the generation of the im- 
ages. (They cover a hexagon centered on the fiber closest 
to the image peak). Figure 2 shows an example of such 
a reconstruction. The visual inspection of these images 
immediately allowed us to identify the loss of lock in the 
guiding system for four exposures, something confirmed 
later by the keywords in the headers of these images. 

The images were adjusted to 2D gaussians with the 
help of the task "n2gaussfit" of IRAF. It is well known 
that the seeing function does not fit a gaussian in the 
outer parts (wings). Therefore, we restrict these fits to a 
relatively small box of 0.8 arcsec (the estimated seeing). 
As a result of these fits we obtained the evolution of 
the following parameters: position (X, Y, position with 
respect to the central fiber), fwhm (semi major axis), 
amplitude, ellipticity, and PA of the major axis. Note 
that for this particular application we are not interested 
in the absolute value of these variables (which depend on 
the assumptions about the function assumed), but only 
on their relative evolution along the night. 

5. RESULTS 

5.1. Auxiliary Variables 

We studied the variation along the night of a set of aux- 
iliary variables, which will be compared with the photo- 
metric index defined in the next section. These variables 
were obtained from the ING meteorological data base, 
from the reduced spectra, and from the reduced images. 
As commented above, an important advantage of the IFS 
(especially the optical fiber approach) is the possibility 
to track independently the variations produced at the 
focal plane of the spectrograph (due, for instance, to in- 
stabilities of the mechanical structures, changes in the 
spectrograph focus, etc.) and those produced at the tele- 
scope focal plane (seeing changes, guiding instabilities, 
sky level variations, etc.). Figures 3-5 show the changes 
of these variables along the night. 

In Figure 3 we show the change in air temperature, 
humidity, atmospheric pressure, sky level, and airmass. 
In this figure we can see that the object was well located 
on the sky for these observations, and we could observe 
it for more than 7 hours at an airmass lower than 2. 
The sky background level was somewhat unstable at the 
beginning of the night. In four images there was a rel- 
atively sudden drop (~ +100 minutes, in the time scale 
of the figure), but it was later discovered that for these 
images we lost the guiding system, so it was a spurious 
result (deviant results from these images were also found 
for other variables — see below). The background has a 
general trend to decrease over the night. The last two 
points in the panel indicate the increase due to the twi- 
light. The values for the air T, humidity, atmospheric 
pressure and wind intensity (not shown) indicate that 
the night had good stable conditions. Measurements for 
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the seeing outside the dome were below 0.8 arcsec most 
of the night. 

Figure 4 shows a set of variables inferred from the spec- 
tra. The two lower panels illustrate the drifts of the 
spectra during the night along and across the disper- 
sion directions (due to the motion of the image of the 
slit with respect to the detector). We can see that the 
spectra move at a relatively constant rate of 0.08 pix/h 
(across) and 0.03 pix/h (along). Although these are very 
small changes, drifts of the spectra over the night make 
the observations sensitive to flat-field corrections, and ex- 
traction aperture location. We did not find appreciable 
relative motion among the spectra, nor global rotation. 
Figure 4 also shows the evolution of the FWHM of the 
spectra integrated over the central band of the index (i.e. 
fi, see below). This also shows a monotonic small change 
during the first part of the night, to reach a stable value 
for the second part of the night. The two deviant points 
(t <~ +100 m.) correspond to the images for which the 
guiding system was lost (the relative change in the spec- 
tra intensity produces a small but appreciable change in 
the FWHM). The variable denoted AFWHM represents 
the relative change in the FWHM between the line (fi) 
and the continuum (/i, / 3 ) bands used for the definition 
of the index as it will be explained below. Although very 
small, we can appreciate a small trend. The variation of 
the total flux integrated over the 37 fibers and over the 
total spectral range used for the definition of the photo- 
metric index (i.e. /1+/2+/3) is also represented, together 
with the individual f\, / 2 , and fz fluxes. As we will see, 
the variation of these fluxes is clearly correlated with the 
star position in the bundle, seeing, and airmass. Finally, 
the top panel represents the ratio of the two continua 
used for the definition of the index (i.e. /1//3). 

Figure 5 shows how the variables obtained after fit- 
ting the reconstructed images with a simple 2D Gaussian 
model change over the night. This figure shows a slight 
drift of the star image over the fiber bundle (~0.4 arc- 
sec) during the first three hours until the point in which 
the auto-guiding was lost, to remain relatively constant 
the rest of the night. Also these panels show a clear 
correlation between the ellipticity and the FWHM, and 
an anticorrclation of these two variables with the ampli- 
tude of the gaussian fit and total flux. This behavior is 
expected due to the fact that, under these good seeing 
conditions, the fiber bundle undersamples the PSF. 



5.2. Index definition 
We define a photometric index probing the NaD lines 

fi - h + h 



as: 



/ = 



fi + h + h ' 



(1) 



where /i,/2, and / 3 represent the sky subtracted flux 
(on a scale of photo-electrons) in the following spec- 
tral ranges: 582.9200-588.5643, 588.5643-590.0887, and 
590.0887-595.7321 nm, respectively. Therefore, the NaD 
feature is well centered within /2, while f\ and / 3 rep- 
resent symmetrical blue and red continua, respectively. 
Note that this index is insensitive to global (i.e. grey) 
variations in the flux. Also note that by increasing the 
width of bands /1 and / 3 (and assuming stable contin- 
uum values) one would reduce the standard deviation 
of a distribution of index values. This apparent gain in 



accuracy is obviously compensated by the fact that the 
index is less sensitive to changes in fi. So, once the width 
of /1 and / 3 are large enough compared with / 2 (so the 
statistical noise fixing the continuum is small compare 
to that due to fi), no real gain is obtained by increasing 
the width of those bands (and further systematic errors 
can be introduced). 

As defined in cq. (1), the index / will probe the NaD 
lines under the assumption that f\ and / 3 are stable con- 
tinua. Charbonncau et al. (2002) have reported very sta- 
ble continuum around the NaD lines, so one could only 
expect temporal variation of f\ and / 3 due to instrumen- 
tal or atmospheric conditions. Some authors (e.g. Winn 
et al., 2004 and references therein) have found relatively 
large time-dependent instrumental (blaze function) vari- 
ations, which were interpreted as caused by flexure of the 
spectrograph. However, in the present case, our fiber-fed 
spectrograph is mounted on a stable platform and, there- 
fore, it is not affected by motions of the rotator. In any 
case, and in order to check instrumental stability, we have 
displayed in Figure 4 the fluxes in the bands fi, fi, and 
/ 3 as a function of time. No apparent differences among 
these three plots can be directly observed, ruling out in- 
strumental variations in flux of the order of 5-10 percent, 
as reported by other authors. In the top panel of Figure 4 
we have also represented the ratio /1//3, which has the 
expected much lower (<1 percent) variation with time 
(i.e. airmass) due the differential atmospheric extinction 
between the two bands. 

From Eq. (1) we can determined how relative changes 
in fi translate into relative changes in /. Specifically, 



AT 

T 



-2(/i + h)f2 



A/ 2 



(/1+/2+/3X/1+/3-/2) h 



and, taking into account that /1 ~ / 3 ~ 4/ 2 , 



Al 



-0.25 



A/2 

h 



(2) 



(3) 



Thus, a relative change of the flux in the NaD band 
translates into a relative change ~ 4 times smaller for the 
index, with the opposite sign. So, taking into account 
that in a typical exposure we collect about 7.15T0 6 e~ 
in fi, the standard deviation due to photon noise for fi 
would be 3.74T0~ 4 and, therefore, for the index / should 
be ~ 9.5T0~ 5 . 

In addition to the index as defined above, we also eval- 
uated a control index using the same definition as in 
Eq. (1), but a relatively line- free band for / 2 . Specifi- 
cally the bands selected for fi, f 2 , and / 3 were: 590.300 
595.146, 595.146-596.454, 596.454-601.300 nm, respec- 
tively. 

5.3. Autocalibration 

The index as defined above was calculated for the 37 
brightest spectra of each exposure. We excluded the four 
exposures for which we lost the the auto-guiding lock, 
and the four exposures taken at airmass larger than 2. 
The resulting 5254 (37 spectra x 142 exposures) values 
are represented in Figure 6 (top) versus the total flux in 
the band (i.e. f t = f\ + f 2 + / 3 ). Clearly there are some 
systematic deviations from the behavior expected in the 
case that noise were only due to photon noise. This is 
particularly obvious for the largest f t values, which are 
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clearly non-linear. At the lowest flux levels the deviations 
cannot be explained by the expected systematic trend in 
a ratio due to statistical errors. They are likely due to 
a small constant background due to scattered light and, 
perhaps, some non-linearity residuals too. In any case, 
these effects can be well calibrated using the data them- 
selves (i.e. autocalibration) . To this aim, we fit all the 
data points to a smooth function (splines3 of order 16) 
of the total flux, Ifu(ft), and we redefine the index as 
I c = I — Ifit(ft)- The corrected photometric index, I c , 
for the 5254 spectra are represented as a function of the 
total flux in Figure 6 (bottom). These new I c values 
should be free of the non-linearity residuals commented 
above, since the dependence with f t has been removed . 

The distribution of values in Figure 6 was also used 
to check and select the proper methodology during the 
reduction process. For instance, we studied the effects 
of resampling the spectra to a common pixel-wavelength 
solution by analyzing the scatter in this figure. It was 
found that the scatter was appreciably larger when the 
spectra were resampled. This plot was also used to check 
the importance of the flat-field and linearity corrections. 
It was found that such corrections do indeed improve 
the quality of the results (i.e. the scatter in Figure 6 is 
reduced), but by relatively modest amounts. 

The next step was to obtain the photometric index 
value for each exposure (time),/'. This was done averag- 
ing the 37 contemporaneous individual values obtained 
in each exposure. We use an inverse variance weighting 
of the index (which is proportional to the total flux, f t ) 
in order to take into account the statistical errors of each 
individual P c value. Therefore, 

r = (4) 

where the subindex i runs over the 37 spectra obtained si- 
multaneously in each exposure, and fl is the correspond- 
ing total flux (in photo-electrons). 

Note that apart from the possibility to perform the 
linearity calibration mentioned above, this procedure for 
generating the index is more advantageous than adding 
all the 37 spectra and generating the index from the com- 
bined spectrum, in terms of the influence of backgrounds 
(i.e. sky, scattered light, etc.) since these are actually 
weighted by the importance of each spectrum for the gen- 
eration of the index at a given time. Also this method is 
also less sensitive to the actual spatial extend used (i.e 
number of fibers involved) . Thus if, for instance, a num- 
ber larger than the 37 fibers considered were used, the 
extra spectra would have had very little weight in as 
a consequence of their low flux. 

Figures 7 and 8 show the dependence of the index with 
some of the auxiliary variables. In Table 1 the degree of 
correlation is indicated in 4th column. Figure 7 and Ta- 
ble 1 show that the strongest correlation (0.76) is found 
against the airmass (which contrasts with the weak an- 
ticorrelation [—0.19], shown by the ' control index', when 
/2 is defined outside the NaD lines, as indicated in Sec- 
tion 5.2). Initially we considered the possibility that this 
strong correlation between the index and the airmass 
could be due to geographical and temporal variations 
of the emission spectrum of the NaD lines produced in 
the Earth's atmosphere. However, the emission spec- 
trum has a very small contribution compared with the 



signal. In fact, considering a typical sky brightness of 
V = 21.8 mag arcsec -1 for La Palma and an EW for the 
telluric emission of NaD of -150 A (B enn and Ellison 
1998; Pedani 2004) , and the parameters used in our ex- 
periment (e.g. width of fi — 15 A, fiber collecting area 
—0.16 arcsec^ 1 , etc.), one should collect more that 10 5 
photons from HD209458 for each photon coming from 
the NaD emission lines. On the other hand, Charbon- 
neau et al. (2002) have shown that the expected drop in 
the NaD band during the transit is ~2.3T0 -4 . So the 
amount we expect to detect during the transit is much 
larger (at least by one order of magnitude) than that due 
to the sky emission (the difficulty of course, is in record- 
ing that signal on top of the photon noise generated by 
the star). Therefore, the strong correlation between the 
index and the airmass is likely to be dominated by the 
absorption produced by the earth atmosphere. The fact 
that the index has a relatively strong correlation with 
airmass (which obviously change with time), may induce 
correlations with other variables which show a temporal 
variation, for which there is not a causal connection. 

5.4. Decorrelation 

By decorrelation we understand the process of remov- 
ing trends in the index correlated with the auxiliary vari- 
ables. This was based on a multiple-regression (least- 
square) fit to the data with a function which is linear 
in the coefficients of independent variables, following the 
prescriptions given by Bevington (1969). Column 4th of 
Table 1 shows the correlation parameter for the 15 se- 
lected auxiliary variables. In order to have an idea of 
the effects of the decorrelation procedure reducing noise, 
we can compare the global standard deviation of the 142 
data values before (1.85- 10 -4 ) and after (1.11 -10~ 4 ) this 
procedure was applied. Note that the standard deviation 
of the decorrelated values is remarkably close to the one 
expected in a purely photon-noise dominated regime. In 
fact, considering a typical level of counts per exposure of 
3.1 • 10 6 ADU (counts) for f 2 , a gain of 2.3 e~/ADU, and 
eq. (3), a standard deviation of the index of —0.95 • 10~ 4 
is expected in a photon-noise dominated regime. (It is 
interesting to note that the fact that the decorrelated 
signal shows a standard deviation close to that due to 
the photon noise suggests that the modal noise is very 
small or it can be well tracked by the selected variables). 

As we can see in Table 1, the airmass is the variable 
which shows the highest correlation with the index. If the 
data are only decorrelated linearly against the airmass 
(i.e. if only the linear trend observed in the bottom panel 
of Fig. 7 is removed) the standard deviation is 1.20 • 
10~ 4 . If we decorrelate against all other variables, but 
we exclude the airmass the standard deviation is 1.17 • 
10~ 4 . In Figure 9 (left) we show the original data (upper 
panel), the fit to a function linear to the airmass (center), 
and the decorrelated data (lower), as a funtion of time. 
Figures 9 (center) and (right) are similar but fitting the 
data to a function linear to 14 (all but the airmass), and 
to all the 15 selected auxiliary variables, respectively. It 
is interesting to note how the rest of the variables also 
carry information about the effects of the airmass. 

If we select for the decorrelation the six auxiliary vari- 
ables derived from the spectra (background, shift x, 
shift y, FWHM, AFWHM, total flux) the standard devi- 
ation is 1.39 • 10 -4 . If we select the nine variables from 
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the images (background, peak, shift x, shift y, distance to 
central fiber (111), FWMH, ellipticity, P. A., total flux) 
the standard deviation drops to 1.24 • 10~ 4 . This illus- 
trates one of the advantages of using IFS: having control 
of the image characteristics at the focal plane offers the 
possibility to improve the photometry by removing cor- 
related noise. 

Although this decorrelation method can be applied 
safely if the functional dependence of the auxiliary vari- 
ables are different to that of the signal to be detected, it 
may destroy part of it when that is not the case. This 
is more likely when many different variables are involved 
in the decorrelation process. Also the fact that our ob- 
servations include only one transit, and this represents a 
relatively large fraction (~30 percent) of the whole ob- 
served period makes this possibility more likely. Ideally, 
observations during long out-of-transit periods will char- 
acterize more accurately the dependence on the variables. 
Obviously having data for multiple transits will minimize 
the possible artificial attenuation of the signal during the 
transit. In any case, to study the degree of destruction 
of the signal during the transit due to the decorrelation 
process we proceed as follows. Similarly to Brown et al. 
(2002) we injected an artificial signal to the data (dur- 
ing the period that the transit took place) before the 
decorrelation was applied, and we analyzed how much of 
this signal was recovered after decorrelation. In a first 
step the effects of each auxiliary variable were analyzed 
individually. That was done comparing the input and 
output signals after decorrelating against each of the 15 
variables. Column 5th in Table 1 indicates the fraction 
of signal recovered, when the decorrelation is done us- 
ing only this variable. In an attempt to estimate the 
error for these conservation factors we used different lev- 
els of injected artificial signal (i.e. 0.5, 1. and 2 times 
the expected change according to the STIS-HST result) , 
but the conservation factors were similar for these three 
cases. Table 1 shows that in several cases more than 
80 percent of the signal is recovered, while in other cases 
this is less than 30 percent. This is the case of three vari- 
ables inferred from the spectra (shift x, shift y, fwhm). 
This is likely due to their monotonic behavior along the 
night, something that can mimic a transit occurred dur- 
ing the first part of the night as it was in the present case. 
However, a variable with a different type of dependence 
such as, for instance, the airmass (which decreases and 
increases along the night, with the inflexion point offset 
with respect to the transit) better preserves the signal. 
During the 132 minutes considered as in-transit observa- 
tions (i.e. excluding the ingress and egress) the airmass 
changed from 1.242 to 1.019. However, after the pass 
through the meridian, the change in airmass from 1.019 
to 1.242 (and on) was done out transit. Then we studied 
the combined effects of several variables. We found that 
when all the 15 variables are used to remove noise (i.e. 
decorrelation process) most of the signal was actually de- 
stroyed: only 17 percent of the signal artificially injected 
is recovered after decorrelation. 

The next step in the analysis was to select a subset 
of auxiliary variables, which lead to a significant reduc- 
tion of noise while preserving a substantial fraction of 
the signal. After several trials based on the individual 
destruction factors, we found that three variables (air- 
mass, FWHM of the image, distance to 111) give a good 



compromise. Specifically, the standard deviation of the 
142 points was 1.19 • 10~ 4 , while 79 percent of the sig- 
nal is preserved after decorrelation. Although the addi- 
tion of FWHM and distance to 111 implies a relatively 
modest improvement in the standard deviation with re- 
spect to the decorrelation against the airmass only, these 
variables were included because they take into account 
the position and shape of the image on the fiber bun- 
dle, which show a relatively strong correlation with the 
total flux (see Fig. 5). In any case, we checked that 
this selection did not affect our results, and values after 
decorrelating against the airmass only will also be given 
hereafter. 

At this point we evaluated the effects of the autocali- 
bration (see Section 5.3) on the final standard deviation. 
If the decorrelation procedure is applied to the data with- 
out applying the autocalibration step a standard devia- 
tion of 1.25 • 10~ 4 is obtained, which compares with the 
value of 1.11 • 10~ 4 mentioned above. (Values for which 
the total flux f t was larger than 1.18- 10 7 , were excluded 
since they show saturation — see Figure 6). This demon- 
strates that, the autocorrelation step allows to improve 
the photometric accuracy as well as to recover some of 
the saturated values (improving the temporal coverage, 
and the error of the mean value). 

5.5. Time Series 

Figure 10 (upper panel) shows the time series data 
values. Vertical lines indicate the contacts. Following 
Charbonneau et al. (2002) we define in-transit observa- 
tions as those that occurred between second and third 
contact (i.e. \t — T c \ < 66.1 m.), and out-of-transit those 
that occurred before first or after fourth contact (i.e. 
\t — T c \ > 92.1 m.). There are 45 in-transit observa- 
tions and 84 out-of-transit observations, ignoring those 
observations that occurred during ingress and egress. 

As mentioned above the standard deviation of all the 
142 values represented in Figure 10 is 1.19 • 10~ 4 . How- 
ever, the standard deviation in-transit (1.46 • 10~ 4 ) is 
significantly larger than out-of-transit (1.03 • 10~ 4 ). We 
think this is attributable to the more unstable condi- 
tions of the atmosphere/background at the beginning of 
the night. Although our meteorological data do not show 
strong evidence for that (see Figure 3), variations of the 
telluric vapor absorption (e.g. Lundstrom et al. 1991) 
and/or of the mesospheric sodium column density (e.g. 
Ge et al. 1998) could produce these instabilities. In any 
case it is remarkable the fact that out-of transit, and con- 
sidering data that expand for a period of ~ 4 hours, the 
standard deviation is very similar to that expected from 
Poison-noise. 

The lower panel in Figure 10 shows the mean in-transit 
and out-of-transit values of the index: (0.12 ±2. 17) • 10~ 5 
and (0.06 ± 1.12) • 10~ 5 , respectively. The quoted l-a er- 
rors correspond to the error in the mean of the statistical 
distribution (i.e. standard deviation/-^). (If the index 
is decorrelated against the airmass only, the correspond- 
ing values are [-0.04±2.16] • 10~ 5 and [0.16±1.13] • 10~ 5 , 
respectively.) These results are compatible with no vari- 
ation of the index during the transit (i.e. the relative 
variation is [0.06 ± 2.44] • 10~ 5 ) and, taking into account 
equation (3), no variation of the relative the flux in the 
NaD lines (i.e. fa) with respect to the continuum dur- 
ing the transit. We also show in this figure the expected 
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change in the index due to a change in ji as the one 
detected by STIS-HST (Charbonneau et al. 2002). Pro- 
vided that these authors found a relative change in f 2 
of (— 2.32 ± 0.57) • 10~ 4 ) during the transit, according to 
equation (3) one should have expected a relative change 
in the index of 5.80 • 10~ 5 with a conservation factor of 1, 
and 4.58 • 10~ 5 with the current conservation factor of 
0.79 (see Section 5.4). (This neglects the slightly differ- 
ent width of the 12 A band used by Charbonneau et al. 
with respect to the 15 A band used by us. This is justified 
by comparing their results for 12 and 38 A bands — see 
Table 2: a change of 3 A in the band width would im- 
ply a change ~ 0.28 • 10~ 5 in the index, about a factor 
of 5 smaller than the error quoted by them). There- 
fore, if a signal with an amplitude as the one reported 
by Charbonneau et al. were real, and after repeating the 
current experiment many times, the (gaussian) distribu- 
tion of results should peak at 4.58 • 10~ 5 with a sigma of 
2.44 • 10~ 5 . However, we (and Charbonneau et al.) have 
performed the experiment only once, and the obtained 
result is about 2 a away from the expected value (for a 
signal as the one detected by Charbonneau et al. data). 
If we take into account the quoted errors, the present 
data and the HST/STIS results disagree by more than 
1 cr, but they are compatible within 1.5 a. 

This limited statistical significance makes us cautious 
to over interpret the present results, which are about a 
factor 2 more uncertain than the HST data. Since at- 
mospheric variables, and especially the airmass, domi- 
nate the decorrelation procedure, a better knowledge of 
the typical time scale for variation of the differential ab- 
sorption in the NaD band by the Earth's atmosphere 
would be useful to further analysis. However, note that 
we could have obtained significantly better accuracy if 
the atmospheric conditions during the transit would have 
been similar to out-of-transit. In that case, the current 
2.44-10 _5 ( 1-er) error would have dropped to 1.9-10 -5 . In 
addition, obviously, multiple transit observations would 
lead also to a significant reduction of this error. The re- 
cent ground based observations by Narita et al. (2005) 
could not confirm nor contradict the Charbonneau et al. 
result. Although these Subaru observations have much 
higher spectral resolution, its lower S/N and temporal 
resolution lead to an error in the predicted relative flux 
in NaD lines (for a 12 A band) larger (by a factor of 
about 4 times) than that of the current observations. 

6. CONCLUSIONS 

The main conclusions of the present work can be sum- 
marized as follows: 

1- We have shown that Integral Field Spectroscopy 
is a powerful technique for obtaining accurate relative 
time series photometry and, therefore, of great poten- 
tial for studying transits of extrasolar planets. The three 
main advantages with respect to previously used meth- 
ods (based on imaging or standard slit spectroscopy) are: 
i) it allows an increase in the duty cycle of the observa- 
tions (and, therefore, the S/N) by distributing the light 
into the two dimensions of the detector; ii) the data can 
be autocalibrated from detector non-linearities and from 
background effects; and iii) since the photometric index 
and the image of the star are extracted from the same 
data-cube, the noise correlated with PSF characteristics 
can be removed. Furthermore, the proposed IFS method 



should not have major demerits with respect to standard 
slit spectroscopy. Although the present study is based on 
ground based observations, the full potential of IFS for 
these type of observations will be reached from space (e.g. 
JWST, TPF). 

2- These observations have also allowed us to explore 
the accuracy limits that can be achieved from the ground 
for these type of observations. In particular, using a pho- 
tometric index that probes the strength of the NaD lines, 
a standard deviation of 1.85 • 10~ 4 was found (during a 
~7 hour period), which is about a factor 2 larger than 
that expected due to photon noise. However, after re- 
moving correlated noise by fitting the data to a function 
linear in three selected auxiliary variables the standard 
deviation drops to 1.19 • 10~ 4 , only ^25 percent larger 
than that due to photon noise. Simulations indicate that 
79 percent of the relative flux change during the tran- 
sit should be preserved after decorrelating against those 
three variables. Our nominal 1-er error for the mean in- 
transit (with respect to the mean out-transit value) for 
the index, is equivalent to 9.7T0~ 5 for the relative flux in 
the NaD band (/2) with respect to the continuum. This 
is substantially better than previous ground observations 
and only a factor ~2 times worse than the 4-transit ob- 
servations with STIS-HST. 

3- The present mean values for the photometric index 
in transit and out-of-transit are compatible with no ex- 
tra depression of the NaD during the transit. Although 
the present results are in apparent contradiction with 
the HST-STIS results, we stress the limited statistical 
significance of this disagreement: they disagree at 1-er, 
but they are consistent within 1.5-cr. We acknowledge 
that our results have larger nominal errors, and are more 
susceptible to systematically underestimate the signal 
through decorrelation. 

4- The present study indicates that a good level of ac- 
curacy can be obtained from the ground with this type of 
IFS observations. The 5254 spectra analyzed here were 
collected during only one night/transit and, therefore, 
are insufficient to reduce the nominal errors to the level 
of the four transit STIS-HST data. Multiple transit ob- 
servations will also allow to analyze the typical time scale 
variations for relevant atmospheric variables and, there- 
fore, better understand the effects of the decorrelation 
procedure on the photometric index. 

5- From an instrument design perspective, this appli- 
cation of IFS can be fostered if: i) the fiber bundle (im- 
age sheer, array of micro-lenses) have an independent 
focussing mechanism, ii) direct illumination of the de- 
tector is possible to generate flat fields, and iii) the sky 
is simultaneously recorded. Fiber systems seem to have 
some generic advantages with respect to another type of 
IFU (e.g. azimuthal scrambling of the light within the 
fibers, the possibility to use static spectrographs, etc.). 
Furthermore, our data indicate that modal noise or the 
temporal variation of the fiber throughput with wave- 
length are not significant sources of noise. 

Further multi-transit observations, with long out-of- 
transit data (ideally from space) should lead to the de- 
velopment of the full potential of the proposed method. 
For ground based observations, the use of the recently 
developed new integral field spectrographs for 8m. class 
telescopes (e.g. GMOS: Allington-Smith et al. 2002, SIN- 
FONI/SPIFFI: Eisenhauer et al. 2003, VIMOS: Le Fevre 
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et al. 2003) also opens an opportunity to push this tech- 
nique forward. 
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TABLE 1 
Decorrelation Parameters 



# 


Variable 


Coefficient 


Error 


Correlation Conservation Factor 


1 


Airmass 


1.074 


5.41E-04 


0.76 


0.88 


2 


Total flux 


8.69E-04 


1.33E-10 


-0.14 


0.88 


3 


Background 


7.74E-02 


-6.61E-07 


-0.48 


0.75 


4 


Spectra Shift (along disp.) 


2.11 


4.64E-04 


0.38 


0.33 


5 


Spectra Shift (across disp.) 


1.60 


3.60E-04 


-0.29 


0.32 


6 


Spectra FWHM 


3.16 


-1.87E-05 


-0.11 


0.38 


7 


Image Total Flux 


5.82E-03 


-1.01E-08 


-0.16 


0.88 


8 


Flux-ratio 


1.68 


-5.12E-04 


-0.32 


0.56 


9 


Gaussian Peak 


4.55E-04 


-6.08E-11 


0.27 


0.77 


10 


Gaussian X 


0.330 


-2.81E-06 


8.9E-02 


0.40 


11 


Gaussian Y 


0.480 


3.20E-05 


-6.3E-02 


0.46 


12 


Distance to 111 


0.384 


3.24E-05 


0.17 


0.99 


13 


Gaussian FWHM 


0.422 


-5.55E-05 


-0.27 


0.84 


14 


Gaussian Ellipticity 


1.59 


6.51E-04 


2.3E-03 


0.85 


15 


Gaussian P.A. 


0.320 


-3.57E-06 


-0.19 


0.94 




All 

1 +12+13 








0.17 
0.79 
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TABLE 2 

NaD measurements of HD209458B 



Work Value Band (A) Experimental setup 



Charbonncau ct al. (2002) (23.2 ± 5.7) X 10" 


-5 


12 


HST/STIS 


(13.1 ± 3.8) x 10- 


-5 


38 




(3.1 ± 3.6) x 10" 


-5 


100 




Narita ct al. (2005) (30 ± 40) X 10~ 


-5 


12 


Subaru/HDS 


Present work (0.2 ± 9.8) X 10" 


-5 


15 
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Fig. 1. — Spectrum generated after combining all the reduced spectra obtained in a single 120 sec exposure. Vertical lines indicate the 
boundaries of the different bands used for the definition of the photometric index (see text). 
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Fig. 2. — Reconstructed image of the star at the focal plane during a typical 120 sec exposure (see text). 
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Fig. 3. — Evolution of some atmospheric variables during the night (see text). Units for the horizontal axis arc in minutes from the 
center of the transit. AU refers to arbitrary units. 
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Fig. 4. — Evolution of some variables inferred from the spectra during the night (see text). From bottom to upper panels: i) Drift of 
the spectra in cross dispersion, ii) shift of the spectra along the spectral direction, iii) FWHM of the spectra, iv) variation of the FWHM 
between the line and continua bands, v) total flux contained in the 37 spectra considered for the definition of the photometric index and 
for the reconstruction of the images, vi), vii, and viii) the corresponding individual fluxes for bands /i, /2, and fy, and ix) ratio /1//3. 
Units for the horizontal axis arc in minutes from the center of the transit. ADU refers to counts on the detector 
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Fig. 5. — Evolution during the night of some variables inferred from gaussian fits to the reconstructed images (see text). From bottom 
to upper panels: i) value of the gaussian peak, ii) distance to the fiber 111, iii) shift along the R.A. direction, iv) shift along the DEC, 
v) FWHM of the gaussian fit, vi) ellipticity of the gaussian fit, vii) P. A. of the major axis, viii) ratio between the flux contained in the 
7 brightest fibers vs. that contained in the 37 selected fibers, and ix) total flux in the images after interpolation. Units of vertical axis in 
ii), iii), and iv) refer to pixels in the interpolated image (1 pixel = 0.05 arcsec). Units for the horizontal axis are in minutes since the center 
of the transit. AU refers to arbitrary units. 




Flux (ADU) 

Fig. 6. — Upper panel. Distribution of the photometric index as a function of the total flux (i.e. ft = /1+/2+/3) for the 5254 spectra. Red 
line indicates the mean systematic behavior. Bottom panel represent the same distribution after removing the mean systematic behavior, 
(see text) 
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Fig. 9. — Temporal variation of the index I c . The upper panels indicate the original distribution of values. Intermediate panels 
represent the fit of the original distribution to a function linear in the variables indicated at top (i.e. left: airmass; center: the 15 variables 
described in Table 1, excluding the airmass; and right: the 15 variables). The lower panels show the decorrelatcd distribution, which is 
obtained subtracting the fit to the original values. The standard deviation for these distributions are: 1.2 • 10~ 4 (lcft); 1.15 ■ 10~ 4 (middlc): 
1.1 ■ 10~ 4 (right). Units for the horizontal axis are in minutes from the center of the transit. 
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Fig. 10. — Upper panel: Temporal variation of the index I c after decorrelation with three selected variables (airmass, image FWHM, 
distance to 111). Vertical lines indicate (from left to right) first, second, third and fourth contact. Horizontal dashed lines show the 
expected ±1-ct according to photon noise. Lower panel shows the averaged values in-transit (i.e. between second and third contact) and 
out-of-transit (i.e. after fourth contact). The value corresponding to the HST observations by Charbonncau et al. (2002) is shown with an 
open square (see text). 



